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Confirming a Population of Hot Dust Dominant, Star 
Forming Ultraluminous Galaxies at High-Redshift 
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ABSTRACT 

We identify eight z > 1 radio sources undetected at 850/xm but robustly detected at 
70/Ltm, confirming that they represent ultraluminous infrared galaxies (ULIRGs) with 
hotter dust temperatures {<Td >=52±10K) than Submillimetre Galaxies (SMGs) at 
similar luminosities and redshifts. These galaxies share many properties with SMGs: 
ultra-violet (UV) spectra consistent with starbursts, high stellar masses and radio 
luminosities. We can attribute their radio emission to star formation since high- 
resolution MERLIN radio maps show extended emission regions (with half light radii 
of 2-3 kpc), which are unlikely to be generated by AGN activity. These observations 
provide the first direct confirmation of hot, dusty ULIRGs which are missed by cur- 
rent submillimetre surveys. They have significant implications for future observations 
from the Herschel Space Observatory and SCUBA2, which will select high redshift 
luminous galaxies with less selection biases. 

Key words: galaxies: evolution - galaxies: starburst - galaxies: high-redshift - galax- 
ies: formation - cosmology: observations 



1 INTRODUCTION 

The most luminous starbursts in the Universe signal a rapid 
growth phase in galaxies. Submillimetre Galaxies (bright 
infrared luminous galaxies with Ssso ^ 2mJy at < z > = 



2.2; e .g. lsinail et al . 2002; Blain et al.ll2002l ; IChapman et all 
I2005D represent some of the best studied ultra-luminous star- 
bursts at high-2, and they are likely building some of the 
most massive galaxies found at the present epoch. How- 
ever, it has been suggested that selecting z —1-3 galaxies 
by their submillimeter emission is prone to fi nding speci- 
mens with cold average dust temperatures (e.g . lEales et al.l 
I2OOOI : iBlain et all |2004 IChapman et al.1 12004 ) leaving the 
possibility of a large population of hotter ULIRGs at high 
redshift. 

The far- infrared (FIR) part of the spectrum (S-lOOO/^m) 
is dominated by the re-processed emission from dust and re- 
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sembles a modified blackbody whose most imp ortant vari- 
able is the characteristic dust temperature, T^ (|Blain et al.l 
|2002| ). At z ~2, the mean redshift of SMGs, submillimetre 
observations at 850/im sample the emission at a strongly 
sloped section of the blackbody, the Rayleigh- Jeans tail, 
where the flux is approximated as Ssso ~ T^^'^ for a given 
Lfir. For example, if Td rises (for fixed Lfir), the peak 
flux shifts towards shorter wavelengths, away from observed 
850^m and the submillimeter (submm) flux can easily drop 
beneath th e submm confusion limit (~2mJy for JCMT; 
iBlainll 19991 ). Current submm instruments restrict studies to 
those ULIRGs with cooler dust temperatures (T^ <40 K) or 
the highest luminosities (L>1O^^L0), potentially missing a 
substantial portion of the luminous starbursts at high red- 
shift. 

Obs ervations of local U LIRGs (e.g. IRAS 60/im star- 
bursts; iLevine et al.l Il998l ) suggest that dust tempera- 
tures can vary substantially in star forming systems, from 
20K up to lOOK. For a IO^Lq local ULIRG, the mean 
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dust temperature is 45K (|Rieke et alj |2009| ) higher than 
36K, the mean temperature of similar luminosity SMGs at 
z ~l-3. A handful of individual ULIRGs at high redshift 
have been found to have characteristic dust temperatures 
around 50-80K (e.g. FSC 10214+4724, Cloverleaf Quasar, 
IRASF15307+3252, APM08297+5255). However, they are 
either exceptionally luminous or strongly lensed, thus de- 
tectable at 70fim, 850fj,m and 1200/im wavelengths with the 
sensitivity of current instruments. This hints that submm 
selection may indeed be biased towards colder dust speci- 
mens, missing some fraction of hotter dust ULIRGs at high 
redshift. 

At 2 ~ 1-3, many galaxies with hot dust could exist but 
are difficult to detect and study as a population. The first 
observa tional effort to charact erize these "hotter ULIRGs" 
came in lChapman et sH (|2004l ). with a selection of submm- 
faint radio-selected galaxies with starburst optical/near-IR 
spectra, similar to SMGs. Since they had similar radio lu- 
minosities, redshifts, and optical spectra to SMGs they were 
thought to be ULIRGs. However, poor constraints on their 
FIR luminosities provided insufficient evidence of their high 
bolometric luminosity. 

We present eight examples of ULIRGs where observa- 
tions suggest hot dominant dust temperatures. ij2] describes 
the data and our observations, |3] discusses the galaxies' 
physical properties and S|4] discusses the selection biases 
prominent at FIR wavelengths and the context of hot-dust 
ULIRGs in other z ^ 2 galaxy populations. Throughout, we 
use AB magnitudes and assum e Ho= 71 kms~^ Mpc^^ and 
^0=0.27 l|Hinshaw et al.ll2009l) . 



2 SAMPLE AND OBSERVATIONS 

We identified 44 potentia l hot high-z; ULIRGs using 
the IChapman et al.l l|2004f ) optically faint radio galaxy 
("OFRG", also submm-faint, star forming radio galaxy, 
"SFRG") selection in the GOODS-N and Lockman Hole 
fie lds. They were detecte d in the ultra-deep VLA radio maps 
of iBiggs fc IvisonI l|2006l ') (Si.4GHz >15/iJy, >3a). Spectro- 
scopic follow up of these potential U LIRGs with Keck LRIS 
revealed starburst spectral features l|Chapman et al ] |2004 
iReddv et al.lliooil '). mostly at redshifts z ^ 1. Obvious spec- 
tral AGN were not included in the sample. There were 26 
GOODS-N sources and 18 Lockman Hole galaxies selected 
by this radio/rest-UV spectra method. Their redshift dis- 
tribution (quartile of z — 1.2-2.6) is similar to the redshift 
distribution of radio-selected SMGs (quartile oi z — 1.2-2.8). 

The Spitzer 70/im maps of GOODS-N and Lockman 
Hole are used for further selection. From the sample of 44 
radio-selected galaxies, five out of 26 in GOODS-N and three 
out of 18 are detected at 70^m above 3a. These eight galax- 
ies make up the sample we discuss in this paper. The posi- 
tional offsets between the 70/^m centroids and the VLA radio 
positions are no more than 5", which is sufficiently smaller 
than the 70fim beam size and astrometric uncertainty. See 
Table [T] for a summary of the observational data. We use 
the remaining 21 GOODS-N 70/im undetected galaxies in 
a stacking test described in section 13.61 (the Lockman sam- 
ple is excluded for this test since its 70^m map has a much 
shallower depth). 

All sources were observed with the Multi-Element Ra- 



dio L i nked Interferometer Network (MERLIN; iThomassonI 
1 19861 : iMuxlow et al] |2005| ), yielding high-resolution 0.3" 
(GOODS-N) and 0.5" (Lockman Hole) beam 1.4-GHz ra- 
dio maps. Three o f the eight galaxies were included in 
the IMuxIow et al] l|2005l ') catalog of HDF radio sources 
(RGJ123711, RGJ123653, and RGJ123645). The 1200^m 
flux limits come from the Max-Planck Millimetre Bolometer 



(MAMB O: iGreve et al. 2008) an d 850Mm flux limits from 
SCUBA (|Borvs et al. 2003 : Coppin et al.|[2006l ). Both fields 
are covered with Spitzer IRAC (3.6, 4.5, 5.8, and 8.0/itm) 
and MIPS (24 and 70^m), however at greater depths in all 
bands in GOODS-N. Optical photometry in GOODS-N is 
from the HST AC^| using the F435W, F606W, F814W, 
and F850LP filters (B, V, i and z bands). The Lockman 
Hole has HST ACS F814W only (PI: Chapman HST 7057). 
One source in Lockman Hole (RGJ105146) is not covered 
by ACS so we use Subaru/Suprime-Cam i-band instead 
(jMivazaki et all 120021 '). The GOODS-N sources have X-ray 
cov erage from the CDF-N. All are formally undetected in 
the [Alexander et al] (|2003l ) catalog above 5a in the full 
band (0.5-8.0 keV), yet four sources are marginally detected 
at fiuxes <6xl0~^®ergs~^ cm~^. Lockman Hole sources 
have X-ray coverage from X M M-Newton (|Brunner et al.l 
|2008| ). and none of the sources are formally detected, 
with a sensitivity limit in the soft band (0.5-2.0 keV) of 
1.9x10"^*^ ergcm"^s"^ and 9x 10"^*^ ergcm"^ s"^ in the 
hard band (2.0-10.0 keV). No X-ray signal is detectable in a 
stack of the Lockman Hole sources. 



3 RESULTS 

3.1 Dust Temperature 

We fit FIR spectral energy distributions (SEDs) to the FIR 
flux densities (at 70/im, 850^m, and 1200/im) using a mod- 
ified blackbody emission curve where 

,,3+/3 



5*1, (X 



(1) 



exp{hv/kTd) - 1 

5*;^ is the observed flux density at rest frequency v,Td is the 
characteristic dust temperature and (5 is the dust's emissiv- 
ity. Since we are limited by a lack of data points in the FIR, 
we must assume that the FIR/radio correlation holds (i.e. 
that LpiR,, evaluated from 8-1000/im, scales linearly with 
1.4G Hz radio luminosity: iCondon.l992i ). While recent works 
(e.g. lBellll2003l : iBeswick et al.ll2008l) indicate that the corre- 
lation might deviate at higher redshifts and lower flux den- 
sities, its effect on the dust temperature calculation would 
be systematic. The 0.2 dex scatter found for the local FIR- 
radio relation suggests an additional uncertainty to our Td 
estimates. We note that the additional uncertainties would 
appl y equally t o the g alaxies in this paper as well as the 
Cha pman et al. I l|2005h SMGs; foremost, we wish to high- 
light the temperature difference between the two popula- 
tions rather than their absolute temperatures. 

^ Based on observations made with the NASA/ESA Hubble 
Space Telescope, and obtained from the Hubble Legacy Archive, 
which is a collaboration between the Space Telescope Science In- 
stitute (STScI/NASA), the Space Telescope European Coordinat- 
ing Facility (STECF/ESA) and the Canadian Astronomy Data 
Centre (CADC/NRC/CSA). 
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Table 1. Observed and derived properties of hot-dust ULIRGs. 



NAME 


z 


Sl.4GHz 


S70 


Sgso 


S1200 


S24 


Lx[2-10kcVl 


LpiR 


SFR 




R-cff 








(MJy) 


(mjy) 


(mJy) (mJy) 


(MJy) 


(ergs-i) 


(Lq) 


(Moyr-i) 


(Mq) 


(kpc) 


(K) 


RG J123710.60+622234.6 


1.522 


38.3+10.1 


3.9+0.5 


<1.8 


<1.8 


227+39 


<9.6xl0"'2 


1.4xl0i2 


247^57 
349^1 


1.2x10" 


1.7 


54+3 


RG J123653. 37+621139.6 


1.275 


86.7+8.3 


6.6+0.4 


<1.4 


<0.7 


164+33 


2.2xl0'*2 


2.0x10^2 


1.7x10" 


2.9 


47+1 


RG J123645. 89+620754.1 


1.433 


83.4+9.8 


4.8+0.4 


<6.3 


<1.6 


172+34 


5.1x10*2 


2.7x10^2 


458tg 


3.0x10" 


2.1 


48+1 


RG J105159. 90+571802.4 


1.047 


74.5+5.6 


7.7+1.2 


<5.9 


<1.4 


738+27 


<2.2xl0''2 


1.0x10^2 


17911* 


6.6x101" 


2.7 


46+2 


RG J105154.19+572414.6 


0.922 


45.4+6.3 


6.1+0.9 


<4.3 


<1.7 


510+22 


<1.6xl0''2 


4.6x10" 




8.3x101" 


2.0 


45+3 


RG J105146. 61+572033.4 


2.383 


33.5+5.8 


9.5+1.1 


<2.4 


<3.1 


298+16 


<1.4xl0"'-'^ 


4.1x10^2 




2.7x10" 


2.1 


72+3 


ULIRGs with AGN: 
























RG J123711.34+621331.0 


1.996 


126.3+8.6 


1.4+0.4 


<6.0 


<4.4 


473+57 


1.9x10*3 


9.7x10^2 




3.5x10" 


3.7 


50+2 


RG J123711 SF component: 




79.6+17.2 












6.1x10^2 


^047+120 

J-'J't' _107 








RG J123649. 66+620738.0 


2.315 


327.0+8.6 


5.8+0.4 


<8.4 


<1.1 


912+78 


1.3x10*^ 


3.7x10^3 




S.lxlQi" 


2.7 


57+1 


B.G J123649 SF component: 




196.6+17.2 












2.2x101=^ 











Observed and derived properties of the eight submm-faint radio galaxies detected at 70^m. SFR and Lpipt are derived from radio 
luminosities, M* from stellar population model fits, and from modified black body FIR continuum fits constrained by 70fim flux 
densities and normalized to radio luminosity. The temperature uncertainties stated here do not take the uncertainty in the FIR/radio 
correlation into account (as described in section [3.1^ . All sources are undetected at 850/im and 1200/im; the fluxes in the corresponding 
columns are 2cr upper limits. Two sources have AGN in addition to luminous star formation (RG J123711 and RG J123649); their LpjR 
and SFR have been adjusted to represent just the star forming components in the abbreviated rows 'RG J... SF component:'; the star 
forming components have been carefully deconvolved from AGN emission in the MERLIN maps as discussed in ^3.41 The effective 
radius, Reff, is the circularized radius corresponding to a surface area of the >3cr MERLIN radio emitting region. 



To calculate temperature, we require an additional con- 
straint to reasonably fit the small number of photometric 
points with a modified blackbody. We fix /3 = 1.5 and adopt 
a single temperature characterization of the emission, al- 
though it is recognized that dust components of a range in 
temperatures are require d to accura tely describe well stud- 
ied nearby galaxies (e.g. iDale et al.ll2001J) . We have tested 
that our fits are relatively insensitive to l3 by also using 
/? = 2, and finding only small (< 2 K) differences in the fit- 
ted Td. Fig. [T] shows the best fit modified blackbodies and 
derived temperatures for our galaxy sample. These black- 
body SEDs are combined with radio synchrotron emission 
(q = 0.8) which is fit through the 1.4 GHz fiux densities. 
Best fit dust temperatures are listed in Table 1. 

3.2 Stellar Mass 

To estimate the galaxies' stellar masses, we combine the pho- 
tometric points in the optical {HST ACS B, V, i, and z for 
GOODS-N and HST ACS i for Lockman Hole) and the 
near-IR {Spitzer-IKAC 3.6/im, 4.5/im, 5.8^m and S.O^m). 
For z =1-3, these photometric points cover the rest-frame 
1.6/xm wavelength where stellar emi ssion peaks. We use the 
HYPERZ photometric redshift code l|Bolzonella et al.l I2OO0I ) 
to fit this photometr y to several stellar population SEDs 
l|Coleman et al.|[l980l ). The optimized stellar population fits 
are shown in Fig. [S] Measured output includes internal ex- 
tinction factor Av and rest-frame K-band magnitude (which 
we use to measure st ellar mass with the method outlined 
by iBorvs et al. I l2005l . assuming a mean light-to-mass ra- 
tio of Lx/M = 3.2). This analysis indicates that the sys- 
tems are quite massive and comparable to the stellar masses 
of the Borys et al. SMGs, wi t h me an ~2xlO"Ma 
We no te that [Alexander et al.l (|2008l ) and IChapman et al.l 
l|2009t 1 point out that the Borys et al. method of using K- 
band luminosities overestimates stellar masses since it does 
not correct for the AGN contribution to 8/im fiux density or 



sample the stellar emission at its peak; while also a probable 
systematic error in our data set, the effect applies equally to 
SMGs and hot-dust ULIRGs. Only one source, RG J123711, 
shows an excess 8/im flux, indicative of AGN power law 
emission; its stellar mass is likely overestimated. Derived 
stellar masses are listed in Table 1. 



3.3 MERLIN Radio Morphology 

We use the high-resolution MERLIN radio maps to assess 
the contribution of AGN to these galaxies by considering 
their radio morphology, which is shown as contour overlays 
on optical imaging in Fig. [3] With resolutions of 0.3" and 
0.5" per beam (for GOODS-N and the Lockman Hole respec- 
tively), the smallest resolvable structure at z~1.5 would be 
3-5 kpc across; however, the typical size of an AGN emis- 
sion region at radio wavelengths is much less than Ikpc. 
This implies that an AGN dominated source w ould be com- 
plete ly unresolved in MERLIN radio maps (e.g. lCasev et al.l 
I2OO9I ). Fig.|3]shows that each of these galaxies has extended 
emission regions on 8 kpc scales, suggestive of spatially dis- 
tributed star formation, unlikely to be generated by AGN. 
The effective radii corresponding to the surface areas of the 
radio emission regions (> 3(t) are given in Table 1. While the 
morphologies are irregular, we find that the effective radii 
average to 2.5±0.6 kpc, which agrees wit h the size measure- 
ments of SMGs in l|Biggs fc IvisonI 120081 ). 

3.4 AGN Content 

As their selection as SFRGs requires, these eight galaxies 
have rest-frame ultraviolet spectral features that are con- 
sistent with starbursts and do not exhibit features of ac- 
tive nuclear activity (e.g. no CIV absorption). Except for 
RG J 123649, their X-ray luminosities (we assume F = 1.8) 
are less than 3x 10*''ergs~i . These low X-ray luminosities 



4 CM. Casey et al. 




Figure 1. FIR to radio spectral energy distributions of eight 
hot-dust ULIRGs. Each galaxy is fit with a modified blackbody 
SED to the 70/im flux density, limits at 850/im and 1200^m, and 
radio luminosity (which is taken to relate directly to Lpia via the 
FIR/radio correlation). On this plot, there is no visible difference 
between /3 = 1.5 and /3 = 2 fits. Inset on each panel are the 
galaxy name, FIR luminosity (integrated from 8-lOOO^m) and 
best fit dust temperature. The two panels marked with "(AGN)" 
show the SEDs for galaxies who show some evidence of containing 
luminous AGN (details described in section [3]4]l . 



Figure 2. Optical and near-IR, photometry points are fit with 
stellar population models. From these SED fits, we interpolate 
rest-K ba nd magnitudes an d compute stellar masses using the 
method of lBorvs et"all ||2005|) : the derived stellar masses are given 
on each panel in the upper left. We also assess the AGN content 
based on an observed 8/^m flux excess, which is only statisti- 
cally significant in the spectrum of RGJ123711. We mark the 
redshifted 1.6/jm stellar bump with a dashed vertical line. The 
two panels marked with "(AGN)" show the SEDs for galaxies 
who show some evidence of containing luminous AGN (details 
described in section [3T4j ■ 



are consistent with star forming Ralaxies with little con- 
tribution from AGN (e.g. lAlexander et al.|[2005h . With the 
exception of RG J123649 and RG J123711, the 6 remaining 
galaxies also have low 24/im flux densities (524 < 250/iJy). 
Low 24/im luminosity is more strongly associated with star 
formation than AGN emission. 

The two hot-dust ULIRGs RG J123711 and RG J123649 
which have the highest X-ray luminosities and mid-IR lumi- 
nosities of the sample might contain powerful AGN despite a 
lack of rest-UV AGN spectral features. The MERLIN radio 
emission in RG J123711 (bottom left panel in Fig.fSjl is sepa- 
rated in two components: an extended region to the NE (63 
percent of the integrated flux), and a bright compact region 
to the SW (37 percent). The X-ray source coincides with this 
SW component. In contrast, a large 10^° Mq molecular gas 
reservoir (a strong case that substantial star formation may 
be present) has been detected in CO and is morphologically 
separated evenly acro ss the two NE and SW components 
l|Chapman et aL I I2OO8I . Casey et al. 2009, in prep). For these 
reasons, we believe the NE region is likely a star forming 



ULIRG merging with the SW region, an AGN-dominated 
galaxy. 

While the existence of an AGN in RG J 123649 is sup- 
ported by a bright X-ray luminosity and a compact core in 
optical imaging, there is no sign of AGN activity in the op- 
tical spectrum or in its ^^8 kpc extended radio morphology. 
We estimate that AGN emission is 40±30 percent and the 
star forming component is 60±30 percent of the radio flux 
by deconvolving an unresolved 0.3" PSF from its radio mor- 
phology. We adjust the radio luminosity and inferred SFRs 
of both RGJ123711 and RG J123649 to solely represent the 
estimated star forming components (see Table 1, "RG J... SF 
component:"). We have also labeled both galaxies as AGN 
in Fig. 1 and Fig. 2 to clarify that they differ from the other 
six. 

It is possible that these two galaxies' 70/im flux densi- 
ties are dominated by emission from torus dust surround- 
ing the AGN (which can have very hot dust temperatures, 
Td ~1000K). However, disentangling the relative contri- 
butions of AGN and star formation to 70/im flux density 
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Figure 3. Contours of MERLIN high-resolution 0.3" radio emis- 
sion overlayed on HST ACS i-band imaging for the eight hot-dust 
ULIRGs. ACS i-band imaging was unavailable for RG J105146 
and we used Subaru j-band in its place. The unresolved MERLIN 
beam size FWHM is shown in the lower right of each panel in 
white. The contour levels plotted are 40, 60, and 80 percent of 
the peak flux density in each image. Structure in radio emission 
is seen in each galaxy on 8kpc scales. 



is beyond the scope of this study since it requires detailed 
knowledge of the FIR SEDs. For this reason, we exclude 
RG J123711 and RG J123649 from further analysis of the hot 
dust ULIRG aggregate properties, since we wish to charac- 
terize the properties of star formation dominated ULIRGs 
exclusively. 

There is no evidence to suggest that the remaining six 
galaxies contain significantly luminous AGN. With a sample 
size of eight galaxies, two of which likely contain AGN, the 



total AGN fraction in the population is 25 percent. This is 
consistent with the estimated AGN fractio n of Submillime- 
tre Galaxies, between 20-40 percent, from [Alexander et al.l 
(,2005. '). 

3.5 Star Formation Rates 

Wc estimate SFRs from the VLA radio luminositi es, using 
the radio/FIR correlation for star forming galaxies (|Condonl 
1 19921 : iHelou et al.lll985l : ISanders fc Mirabellll996l ). We ab- 
stain from calculating UV-inferred SFRs since they would 
be subject to significant (yet uncertain) extinction factors 
(as evidenced by the irregular MERLIN radio morphologies 
relative to the distribution of rest-UV fiux in Fig. [Sjl . The 
FIR luminosities derived from the radio are consistent with 
the submm detection limits and the fiux densities at 70/im as 
shown in Fig. [T] The inner quartile (25-75 percent) of star 
formation rates in this sample is 200-1100Mq yr~^ (both 
derived quantities, Lfir and SFR are given in Table 1). 

Star formation rate densities are found by dividing these 
SFRs by the surface areas of MERLIN radio emission re- 
gions. Comparing these SFR densities to their theoretical 
maximum— the maximum gas den sity divided by t he local 
dynamical time (see equation 5 of lElmegreenI [l999l . we use 
idyn~ 4xl0^yr)— we can determine if the implied SFR den- 
sity exceeds the theoretical prediction. While local ULIRGs 
with EspR «200Mq yr~^ kpc~^ are forming stars at their 
theoretical maximum (e.g. iTacconi et al]|2006l ). none of the 
galaxies come within a factor of four of exceeding their max- 
imum SFR density limits. 

We note that one of the eight galaxies of our sample, 
RG J1 23710, has been detected in CO gas by iDaddi et al.l 
(|2008t ). The presence of a molecular gas reservoir >10^° M© 
indicates that high levels of star formation (>100Mq yr~^) 
can potentially occur. In their interpretation of the detec- 
tion, Daddi et al. claims that RG J123710 is most likely a 
spiral galaxy undergoing modest-efficiency star formation 
lower than most SMG ULIRGs. However, our detection of 
RG J123710 at 70^im directly confirms that it is a ULIRG. 
More CO observations are needed to investigate if the rest of 
the hot-dust ULIRG population exhibits modest star forma- 
tion efficiencies or high, SMG-like efficiencies (Casey et al., 
in preparation). 

3.6 70^m Undetected SFRGs 

Our eight galaxy sample likely represents the most lumi- 
nous subsample of the hot-dust ULIRG population because 
of the current limitations in 70/im surveys. To investigate 
the possible lower-luminosity extension of the population, 
we analyze the 21 GOODS-N galaxies which were selected 
as SFRGs but not detected at 70/im (the Lockman sample 
is excluded since its 70^m coverage has a much shallower 
depth). Here, we stack 70^m cutouts of these 21 SFRGs, cen- 
tered on the ir VLA positions. W e us e the stacking method 
described in iHuvnh et al.l (|2007l ) and iHainline et al.l (|2009l ) 
which was used to stack ~68 70^im undetected SMGs. We 
have a detection of 0.45±0.15 mjy in the SFRG stacked im- 
age at 3.0(7, which is comparable to the 0.48±0.15mJy 3.2a 
detection for the SMGs of Hainline et al. This indicates 
that SMGs and 70/jm-undetected SFRGs likely have sim- 
ilar 70/im luminosities, although larger samples of SFRGs 
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are needed to see if their relative detection rates at 70/im 
are similar. Because of their detection at 70/im, the eight 
SFRGs of this paper are likely the highest luminosity (z<3) 
and potentially the hottest dust temperature specimens of 
the SFRG population. 



4 DISCUSSION 
4.1 Volume Density 

The volume density is estimated by constraining the density 
of galaxies satisfying our selection criteria. The density of 
these 70fim detected ULIRGs is 1.4x10"^^ Mpc"^. However, 
the difference in volume density estimates between GOODS- 
N and Lockman Hole, which is expected due to the 70/im 
sensitivity differences between the two fields, prevents us 
from constraining the hot-dust ULIRG luminosity function 
well. 

While this explicit population of hot-dust ULIRGs is 
~ 5x more rare than SMGs, there likely exists a larger 
population of z > 1 ULIRGs with hot dust temperatures. 
Firstly, spectroscopic incompleteness of ~50 percent sug- 
gests the sample is likely twice as large (although it is noted 
that SMGs are also affected by spectroscopic incomplete- 
ness). Secondly, there are likely additional ULIRGs a.t z ^ 2 
undetected in the radio (for instance from scatter in the 
FIR/radio relation). Finally, the 70/xm sensitivity limits of 
these maps do not exclude many of the other SFRGs in 
our parent sample from having warm enough dust to still 
represent hotter ULIRGs than SMGs. As we have shown in 
our analysis of the GOODS-N sample, proportionally more 
SFRGs are 70/im detected than SMGs; as also indicated by 
the stacked result of 70pim undetected SFRGs, many more 
SFRGs might become 70Atm detected if the GOODS-N 70^m 
depth can be achieved in other fields. The Herschel Space 
Observatory will dramatically improve the uniformity and 
depth of 70^m coverage (confusion limit of 0.04 mjy at 2a), 
potentially discovering many of the ULIRGs that lay just 
below the current MIPS detection limits. 



4.2 Selection Effect in Lpm-Td and Redshift 

Fig. |4] illustrates the selection effects of IR- wavelength ob- 
servations as a function of FIR luminosity and characteristic 
dust temperature. Assuming the radio/FIR correlation, ra- 
dio detection is independent of dust temperature. At z = 2, 
a 10^^ Lq system is only detectable at 850/im or 1200/im if 
it has Td < 40 -ff. Since the mean redshift of radio-detected 
SMGs is 2; = 2.2, the dust temperature selection-effect sig- 
nificantly impacts the proposed surface density of ULIRGs 
and thus their contribution to the global star formation rate 
density near 2 ~ 2. If instead we select objects by the meth- 
ods of this paper at 70/im, the dust temperature selection 
effect is inverted, and detection is possible for the hottest 
z ~2 sources (T^ > 7GK). By lowering the redshift bound- 
ary from 2 > 2 to z > 1, detection is possible for 10^'' L0 
galaxies that have T^ > 35 if. 

While the FIR selection functions have hard boundaries 
in LpiR and Td, their dependence on redshift is important. 
From Fig. |4l we see that 70//m selection is not ideal for 
z > 2 samples where it is limited to temperatures Td ^ 70 K 
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Figure 4. Far-IR luminosity against dust temperature illustrates 
the temperature/luminosity-biases of different single-wavelength 
detection techniques, at 70/im (red), 350/im (magenta), 850/im 
(blue), 1200/tm (green), and 20cm/1.4GHz radio (gray). The 
boundaries on the plot are determined by flux density 2ct upper 
limits for these observations: 1 mJy at 70 /xm, lOmJy at 350 /xm, 
Smjy at 850 /^m, 1.5 mJy at 1200 /im, and 10/tJy at 20 cm (the 
limits for GOODS-N). The boundaries also represent a limit in 
redshift which is chosen based on the characteristic mean red- 
shift of sources detected at that wavelength. Sources may be de- 
tected at the given wavelength on the high-luminosity side of each 
boundary. We overplot the SMG population which is selected 
at 850/xm l|Chapman et al.l I2OO5I : iKovacs et aP 12006^1. and z~2 
1200/im-luminous ULIRGs from the EGS field llYounger et al] 
i 200d ). and the sample in this paper, detected at 70/im with a 
2 = 1 boundary (solid) and 2 = 2 boundary (dashed). 

whereas at z ~ 1 it can select a wider range of ULIRGs with 
Td ~ 40 K. In Fig. Owe show the projections of both Lfir 
and Td with redshift. This highlights that the luminosities 
of the 70/im-selected sample are comparable to SMGs at the 
same redshift. However, as a function of dust temperature, 
the 70/im galaxies occupy an unexplored region in z — Td 
which is heavily restricted by the selection biases illustrated 
in Fig. m 



4.3 Comparison of BED shape to SMGs 

A direct comparison of the observationally derived SED 
for these hot dust ULIRGs with SMGs is shown in Fig. |6l 
We average the 70/im detections, 850/im upper limits, and 
radio detections of the six star formation dominant hot 
dust ULIRGs to form our composite SED; the composite 
SED adopts the mean redshift {z = 1.5), FIR luminos- 
ity (1.9x10^^ L0), and dust temperature (52 K) of the six 
galaxy sample. Since SMGs have a mean redshift of 2.2, 
we take several comparable sub-samples of GOODS-N and 
Lockman Hole SMGs to generate an average SMG SED 
with mean redshift 1.5. The resulting mean luminosity of 
this SMG sub-sample is 2.2x10^^ L0 with dust temperature 
36 K. 

Based on these SED fits, we do not expect SMGs to be 
detected at 70/im. This agrees with observations; nearly all 
SMGs (68/73, >93 percent) are undetected <3a at 70/im 




0.1 1.0 

REDSHIFT 

Figure 5. Far-IR luminosity and dust temperature as functions 
of redshift (log scale), with the same symbols as in Fig.|4] with a 
low redshift 350/^m galaxy population overplotted for comparison 
dYang et al.ll200^ . Above, the 70/im has similar luminosities to 
SMGs at the same redshifts. Below we highlight the poorly ex- 
plored region (shaded region) in dust temperature with redshift. 
The 70/im-selected sample begins to sample this unexplored area. 
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Figure 6. A composite spectrum of the six star forming, 70/im 
detected, hot dust ULIRGs relative to a composite SED for a 
sub-sample of SMGs with < z >=1.5. We select several sub- 
samples of SMGs all with mean redshift 1.5 to mimic the redshift 
distribution of the hot-dust ULIRGs. The average flux densities 
for hot-dust ULIRGs are 6.4 mjy at 70/im, <2mjy at 850/im and 
50/iJy at 20 cm. The flux density points for SMGs are <1.5mjy 
at 70/im, 7mjy at 850/tm and 60/tJy at 20cm. The red and blue 
vertical bands illustrate the range of peak fluxes corresponding 
to two temperature regimes; blue represents 36ib7K and submm 
detected, while red represents 52±10K. Note that at 350/tm, both 
samples of 70/im and 850/tm galaxies would be equally easy to 
detect. 



iHainline et al.ll2009l ). Those that are 70/im-detected are at 
low- 2; and have higher luminosities than most SMGs at low 
redshifts. To compare the SMG and SFRG populations, we 
compare their 70/im detection rates. We limit this test to 
z < 2 since the redshift selection functions for the popula- 
tions differ (submm-|-radio s election and radio sele ction have 
different redshift biases; see IChapman et al.ll2005l ). This re- 
moves most effects of redshift-bias from radio detection in 
both samples. SMGs and SFRGs have similar radio lumi- 
nosities, so the same fraction of each should be detected at 
70/im if they have the same distributions in dust tempera- 
ture. In GOODS-N, where the observations are deepest and 
the most complete, 4/8 (50±25 percent) z<2 submm-faint 
radio galaxies are detected at 70/im but only 1/14 (7±7 
percent ) of submm-bright radio galaxies (SMGs) a re de- 
tected ((chapman et all (20051 ; fPope et al.11200^ . |2008| ). The 
Lockman Hole observations do not have nearly the depth as 
those of GOODS-N and thus have fewer sources; however, 
the statistics are consistent with the finding in GOODS-N, 
with 3/7 (43±25 percent) submm-faint galaxies detected at 
70/im and only 1/5 (20±20 percent) SMGs detected. There- 
fore, in contrast to SMGs, galaxies selected by our method 
(SFRG selection) are far more likely to be detected at 70/im 
(30-50 percent). 

4.4 Comparison to other IR Selection Techniques 

Since the advent of Spitzer and other infrared observatories, 
the population of dusty galaxies studied at high redshift has 



grown substantially. Many of these galaxies are more bolo- 
metrically luminous than typical UV-selected galaxies, but 
they may be poorly understood as a population due to faint- 
ness in the optical. Furthermore, the optical counterparts 
for 850/im and 70/im sources are not easily identified due to 
the ~12-20" beam size; for these reasons, substantial por- 
tions of SMGs and ULIRGs have no redshift identification. 
Similar problems exist for other IR-selected galaxies. Here, 
we contrast the hot dust ULIRGs in this paper to other 
infrared-selected, dusty galaxy populations at 2 ~ 2, noting 
potential overlap of selection techniques. 

To circumvent the problems that arise from FIR selec- 
tion (increased beam size and poor sensitivity limits) , many 
dusty galaxies are selected in the mid-IR by their observed 
24/im flux densities. However, the selection of dusty star- 
forming LIRGs and ULIRGs through mid-IR 24/im contin- 
uum diagnostics is frequently contaminated by power law 
emission from AGN, and inferring star formation properties 
can require large correction factors to bolometric luminosity 
(e.g. Dale ct al. 2005; Papovich ct al. 2005). All eight galax- 
ies in o ur sample are classified as Dust O bscured Galaxies 
(DOGs; [Dev et al.ll2008l : |Pope et al.ll2008l ): however, the se- 
lection criteria of DOGs is so broad that without detailed 
additional information (the radio maps, FIR detections and 
non-detections, and near-IR photometry) it is difficult to un- 
derstand the importance or evolutionary significance of this 
classification. We further emphasize that 3/8 specimens in 
our 70/im sample have S24 <300/iJy (and are only selected 
as DOGs under the selection described by Pope et al.), re- 
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vealing that ultraluminous activity at z =1-3 can be missed 
by 24/im selection criteria; 

A recent study by lYounger et alJ (|2009l ) presents 12 
2 ~ 2 ULIRGs in the EGS field, selected in the Spitzer 
IRAC and MIPS bands and foUowed up with MAMBO at 
1200^m (where 9/12 are detected). They have a character- 
istic dust temperature 41 ± 5 if (using /3 = 1.5 modified 
blackbody fits) and are shown in Fig. |4] and Fig. [5] The 
12 Younger et al. galaxies bridge the gap between SMGs 
and hot 70/im-detected sources in temperature but are much 
more luminous than the 70/im sample. They have a mean 
1200pm flux density of ~ 1.6 mjy where our 70/im sam- 
ple is uniformly undetectetfl. This would imply that the 12 
EGS galaxies have 850/im flux densities ~3-6mJy, FIR lu- 
minosities (8-1000/im) >1O^''L0 and thus could be classi- 
fled as SMGs. The selection of the Younger et al. sample in 
the mid-IR, by the presence of a prominent stellar bump in 
IRAC photometry and dust at 24/^m, is similar to the DOG 
selection, although with a few additional constraints reveals 
a potentially strong probe of ultraluminous activity at high 
redshift. However, more far-IR observations of z ~ 2 mid-IR 
selected samples are needed to assess the selection's com- 
pleteness in choosing both star forming ULIRGs and AGN 
dominated ULIRGs at z > 1. 

Selection at 350pm at z~2 has similar temperature 
biases as 850pm and 1200pm due to its high flux limit 
(S350 >10mJy at 2a). As seen in Fig. |S1 350pm selection 
at z ~ 1.5 has the potential to select SMGs and hot-dust 
ULIRGs equally well since it samples the dust SED near 
its peak for a wide range of temperatures. However, until 
the depth of 350pm imaging can be improved, current ob- 
serving facilities allow only 350pm dete ctions of the mos t 
luminous or low redshift sources (e.g. lYang et ah I I2OO7I ). 
The Balloon-borne Large Aperture Submillimeter Telescope 
(BLAST) has recently mapped the Chandra Deep Field 
South at 250pm, 350pm, and 500pm at greater depth than 
has been done before in any fleld at those wavelengths 
(|Patanchon et al]|2009l ) . While useful for measuring the flux 
densities of already known ULIRGs and high-redshift galax- 
ies, the statistically significant source counts in the BLAST 
map are limited (with only ~12 galaxies detected at z>l; 
IViero et al.ll2009D . When it begins full operations, the Her- 
schel Space Observatory will be another observing tool at 
350pm and will dramatically improve the depth of 70pm 
observations in many heavily observed fields, thus providing 
better statistics at shorter FIR wavelengths. In addition, fu- 
ture work from the SCUBA2 instrument at 450pm have the 
potential to expand the sample of 2: ~ 1 — 3 ULIRGs through 
deeper observations near the peak of the dusts' SED. 



5 CONCLUSIONS 

This paper has observationally demonstrated that submil- 
limetre wavelength observations of galaxies at z ~ 1 — 3 have 
a strong temperature bias which preferentially detects galax- 
ies with cooler dust temperatures. By selecting ULIRGs that 
are detected at 70pm (and verifying that they are starburst 

2 The depth of MAMBO coverage in GOODS-N is ~0.7-0.8mJy 
RMS compared to the ~0.4mjy photometric RMS of the EGS 
pointings, while Lockman Hole has 0.8-1.0 mJy RMS. 



dominated) we have found a set of galaxies which are un- 
detected in the submillimetre yet still represent some of the 
most luminous systems at 2: > 1. Their volume density is 
1.4x10"^ Mpc~^, approximately 5x more rare than SMGs, 
although this has likely been underestimated due to spectro- 
scopic incompleteness and large variations in the sensitivity 
limits of 70pm surveys. 

Like most IR-selected dusty galaxies, the SFRs in these 
hot dust ULIRGs are on average much higher than the 
majority of UV-selected star-forming galaxies at the same 
epoch. While some UV-sel ected sources have high stellar 
masses (jShaplev et al ] |2005l , Hainline et al., in preparation) 
they show consistently lower SFRs, reiterating that the UV 
misses the most dramatic star-forming galaxies at high red- 
shift. Due to their similarity with SMGs (in terms of stellar 
mass, UV-spectra, FIR and radio luminosities), the hot-dust 
ULIRGs represent an extension of SMGs towards a wider 
range of dust temperatures at z~ 1—3, similar to the al- 
ready observed wide range in dust temperatures seen in local 
ULIRGs. 

We have shown that AGN contribute little to the FIR 
luminosities in six of the eight hot dust ULIRGs. The de- 
tection of extended radio emission (with typical radii 2- 
3kpc) suggests spatially distributed star formation rather 
than compact AGN. In addition, the UV-spectral properties 
and low X-ray luminosities are consistent with star forma- 
tion. The two ULIRGs that do show signs of AGN activity 
also have evidence for substantial star formation and are 
therefore comparable to SMGs which contain bright AGN. 

Combining our hot-dust ULIRG population with the 
SMGs results in a sample of z ~ 2 ULIRGs with less 
temperature-dependent bias than the SMG population 
alone. However, we have shown that searches for high red- 
shift FIR luminous galaxies are likely still somewhat in- 
complete. Even at 2 ~ 2, there are likely other FIR lu- 
minous galaxies which are not well characterized by cur- 
rent submm or FIR observations. Herschel and SCUBA2 
will help push detection of more z > 1 ULIRGs in the 50- 
500pm range while more work with deeper radio observa- 
tions (from e-MERLIN and EVLA) will be needed to find 
potential ULIRG activity beyond z > 3. By working towards 
completeness in the z > 1 ULIRG population, we will learn 
about the role of heavy, short-lived star formation in the 
formation and evolution of galaxies in the early Universe. 
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